disp('gamma distribution with parameter alpha, beta')
alpha = input('parameter alpha:');
beta = input('parameter beta:');
x = linspace(0,10);
plot(x, (beta^alpha / gamma(alpha)) .* x.^(alpha-1) .* exp(-beta.*x));
title(sprintf('gamma distribution with alpha=%0.5f beta=%0.5f', alpha, beta));
graph_name = sprintf('gamma_%0.5f,%0.5f', alpha, beta);
print(graph_name, '-dpng');
close(1);
